Nucleon Instability in SUSY Models 



PRAN NATH 

Department of Physics, Northeastern University, Boston, 
MA 02115, USA 

R. ARNOWITT 

Center for Theoretical Physics, Department of Physics 
Texas A & M University, College Station, TX77843, USA 

A review of the current status of nucleon stability in SUSY unified models is given. 
The review includes analysis of nucleon stability in the minimal SU(5) supergravity 
model, and in extensions of the minimal model such as SO(10) models, and models 
including textures. Implications of the simultaneous constraint of the experimental 
limit on the proton lifetime and the constraint that dark matter not overdose the 
universe are also discussed. 



1. Introduction 

This review concerns the status of nucleon stability in supersymmetric 
theories. First, we will give a general discussion of baryon number violation 
relevant to proton stability in SUSY theories. Next we will discuss in some 
detail the predictions of nucleon stability in the minimal supergravity model. 
Finally we will discuss proton stability in extensions of the minimal model 
and comment on the issue of proton stability in string models. There are, of 
course, many sources of baryon number violation in unified theories. First, 
grand unifed theories have lepto-quark mediated proton decay and in SU(5) 
the proton decay via lepto-quarks is given bja 

T(p^e+TT°)^( ) 4 10 31± V (1) 

In SUSY SU(j5) models one estimates My = 1.1 x 10 16 GeV, which gives a p 
decay lifetimtn of 1 x 10 35±1 . The current limit experimentally ifl 

r(p -> e + 7T°) > 9 x 10 3 V, (90%CL) (2) 

and it is expected that Super KamionkandelSuper-K) and Icarus will be able 
to probe the e + 7r° proton decay mode up tcfl 1 x 10 34 y. Thus the e + 7r° mode 
in SUSY SU(5) may be on the edge of detection if Super-K and Icarus reach 
their their maximum sensitivity. 

In supersymmetric theories there are sources of proton instability arising 
from terms in the superpotential. The main purpose of this talk is to review 
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the status of proton stability in the presence of these additional sources. The 
outline of the rest of the paper is as follows: In Sec. 2 we discuss B and L 
violation in supersymmetric theories from the superpotential terms and show 
that generically even with R parity invariance there is p decay from dimension 
five operators in SUSY/string unified models. In Sec. 3 we discuss p decay in 
the minimal SU(5) model. In sec. 4 we discuss p decay in extensions of the 
minimal model. Conclusions are given in Sec. 5. 

2. Sources of P Decay in SUSY/String Models 

We turn now to a discussion of the sources B and L violating interactions 
in the superpotential and their effect on proton stability. It is well known 
that B and L violating dim 4 operators lead to fast proton decay. Thus, 
for example, the effective dimension four low energy SU(3) x SU(2) x U(l) 
invariant interaction given by 

W = X u Qu c H 2 + XdQd ^ + X e Le c H x + m#i#2 

+ {X' B u c d c d c + X , L Qd c L + X'[LLe c ) (3) 

The terms in the bracket lead to fast proton decay and consistency with ex- 
periment requires the constraint 

X'b X'l < O(10- 27 ) (4) 

This type of decay is eliminated in the MSSM by the imposition of R parity 
invariance [where R — (_i) 3B + L + 2S ]. it is likely that this discrete R symmetry 
is remnant of a global continuos R symmetry. In that case there is the danger 
that the global symmetry may not be preserved by gravitational interactions. 
For example, worm holes can generate baryon number violating dimension 4 
operators and catalize proton decajd. A decay of this type is not suppressed 
by either a large exponential suppression or by a large heavy mass and is thus 
very rapid and dangerous. 

One way to protect against wormhole induced fast proton decay of the 
type discussed above is to elevate the relevant global symmetries which kill 
B and L violating dimension 4 operators to gauge symmetriefHj. In this case 
if the local symmetry breaks down to a discrete symmetry, then the left over 
discrete symmetry will be sufficient still to protect againt worm hole type 
induced proton decay. There is, however, the problem that in unified theories 
undesirable proton decay can arise as a consequence of spontaneous symmetry 
breaking even if the dangerous dimension 4 operators were forbidden initially. 
Thus, for example, for SO(10) one can have a interaction of the type (16) 4 
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which gives terms of the following type in the superpotential 



u c d c d c v c , QLd c u c (5) 

mpi mpi 

If there is a spontaneous generation of VEV for the v c field, then terms of the 
type u c d c d c and terms of the type QLd c emerge which once again lead to a 
rapid proton decay unless < v c > /Mpi is 0(1CU 13 ). Thus one must make 
certain that higher dimensional operators in the superpotential do not lead to 
dangerous p decay after spontaneous symmetry breaking takes place or that 
spontaneous symmetry breaking does not occur. 

Next we discuss p decay from dimension 5 operators (dimension 4 in the 
superpotential) which contain B and L violating interactions. Using SU(3) x 
SU(2) x U(l) invariance one can write many B and L violating interactions. 
Examples of such terms in the superpotential involving matter fields are 

QQQL, u c u c d c e c (6) 
One can also have terms which involve the Higgs, e.g., 

QQQH ll Qu c e c H 1 ,etc. (7) 

The second set of terms can be eliminated if we impose R parity invariance. 
However, B and L violating terms of the first type arise quite naturally in 
SUSY unified models via the exchange of color Higgs triplet fields. In fact 
we now show that most SUSY/string models will exhibit p instability via the 
Higgs triplet couplings. Tha,p decay arising from the color Higgs exchange is 
governed by the interactional 

HiJ + KHi + HiMijHj (8) 

where J and K are quadratic in the matter fields. It is then easily seen that 
the p decay from dimension five operators is suppressed provided 

(M-% = (9) 

Now a condition of this type can be met in one of the following two ways: 
(i) discrete symmetries, and (ii) non-standard embeddings. However, most 
SUSY/string models do not normally have discrete symmetries of the desired 
type which automatically satisfy Eq.(9), and only very specific SUSY/string 
models (the flipped models!!) satisfy (ii). Thus in general SUSY/string models 
do not have a natural suppression of p decay via dimension 5 operators. Thus 
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suppression of p decay here must be forced by making the Higgs triplets heavy. 

3. Proton Decay in Minimal SU(5) SUGRA Model 

We begin with a discussion of the comparison of the minimal SU(5) model 
predictions for the gauge coupling constants at the electro- weak scale with the 
LEP data. It is well known that the model with the MSSM spectrum extrapo- 
lated to high, energy gives a reasonably good fit for the coupling constants with 
experimental. However, a more accurate analysis shows that the_thcoretical 
value of a s predicted by SU(5) is about 2a higher than experimenfiilO. There 
are a variety of ways in which one can achieve a correction of size 2a. These 
include Planck scale corrections, and various other extensions of the minimal 
SU(5). We discuss here briefly the possibility of Planck scale corrections which 
are expected to be of size 0(M/ Mpi) where M is size of the GUT scale. It is 
reasonable to expect that such corrections are present due to the proximity of 
the GUT scale to the Planck scale. Effects of this size can arise via corrections 
to the gauge kinetic energy functiorE30'L3. For example, for the case of SU(5) 
one can introduce a field dependence in the the gauge kinetic energy function 
scaled by the Planck mass so that 

f a p = (5af3 + ^j^ da ^ Yjl ) ( 10 ) 

Here c parametrizes Planck physics and E is the 24-plet of SU(5). The analysis 
shows that it is easy to generate, a 2a correction to a s with a c ~ 1 to achieve 
full agreement with experiment! It is also possible to understand a 2 a effect 
on a s from corrections arising from extensions of the minimal SU(5) modeLsuch 
as, for example, in some modified versions of the missing doublet modeiJ. 

Computation of the Higgsino mediated proton decay lifetime in super- 
symmetric theories involves both GUT physics as well as physics in the low 
energy region via dressing loop diagrams which convert dimension five oper- 
ators into dimension six operators which can be used for the computation of 
proton decay amplitude£3'L3. Now the dressings involve 28 separate sparticle 
masses and many trilinear couplings which mix left and right squark fields. 
Thus is general no quantitative predictions of proton decay can be made in the 
MSSM which has a large number of arbitrary parameters in it. In the minimal 
SUGRA modetJo the number of arbitrary parameters is vastly reduced. Us- 
ing the constraints of radiative breaking of the electro- weak symmetry one has 
only four arbitrary parameters and one sign in the SUSY sector of the theory. 
These can be chosen to be 
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too, mi , Aq, tan/3, sign(u) (11) 

where too is the universal scalar mass, mi, is the universal gaugino mass, Aq 
is the universal trilinear coupling, tan (3 =< H 2 > / < H\ >. Here Hi gives 
mass to the down quark and H 2 gives mass to the up quark, and u is the Higgs 
mixing parameter. Because of the small number of parameters, the minimal 
supergravity model is very predictive. In turns out that as a consequence of 
radiative breaking of the electro-weak symmetry, over most of the parameter 
space oLthe minimal supergravity model one finds that scaling laws hold and 
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ha£3 



~mg (fi < 0);to^ - jm § {a > 0) 
2m 2l ~ to^ ~ to^; to^ 3 ~ to^ 4 ~ m^ 2 » m Si 

m° H ~ niA ~ to^± >> to/j (12) 

In the following we shall use the framework of supergravity to compute 
proton decay amplitudes. For concreteness we shall use the minimal SU(5) as 
the GUT group and extensions to the non-minimal case will be discussed in 
the next section. 

The interactions which govern proton decay in minimal SU(5) are given 

by 

W Y = -^f uj e uvwxy H?My w M* v + f 2lJ H 2u M lv M™ (13) 

After breakdown of the GUT symmetry and integration over the Higgs triplet 
field the effective dimension five interaction below the GUT scale which governs 
p decay is given bjJi 

L L 5 = -^e abc (Pf? V) l3 (f 2 d ) kl (ii LM d Lcj (e c Lk (Vu L ) al - v c k d Lal ) + ...) + H.c. 

= --^e abc ( V^r)^{PVf d ) kl {e c Rl u Ra] u Rck d Rbl + ...) + H.c. (14) 

where Lg is the LLLL dimension five operator and is the RRRR dimension 
five operator. The Yukawa couplings can be related to the quark masses at 
low energy by 



= ft (sin29 w /e)M z sin/3 
mf = f?{sin26 w /e)M z sinf3 (15) 



5 



and Pi are the inter generational phases given by 

Pi = (e*»), $^7i = 0; »=1,2,3 (16) 

i 

There are many possible decay modes of the proton. The most dominant of 
these are those which involve pseudo-scalar bosons and leptons. These are 

ViK + ,Vi-K + ;i — e,fi,T 
e+K°, n+K°, e+7r°, fi+ir°, e+r), (17) 

One can get an estimate of their relative strengths by the quark mass factors 
and by the CKM factors that appear in their decay amplitudes. These are 
listed in Table 1. Additionally the decay amplitudes are governed by relative 
contributions of the third generation vs second generation squark and slepton 
exchange in the loops. The relative contribution of the third generation vs 
second generation exchange in the loops is governed by the ratios y tK , y tv ,..etc. 
These are also listed in Table 1. 



Table 1: lepton + pscudoscalar decay modes of the proton 



SUSY Mode 


quark factors 


CKM factors 


3rd generation enhancement 


v e K 


mdrric 


V? 1 V 2 iV 22 


(l + yl K ) 


v„K 


m s m c 


v 2 \v 21 v 22 


(i + y l 2 K ) 


v T K 


m b m c 


V^V 21 V 22 


(i + y$*) 




m d m c 


vl x v 2 \ 


(l + yf) 




m s m c 


v 21 v 21 


(i + yf) 


u t tt, v r r\ 


m b m c 


v 31 V 21 


{i + yf) 


eK 


m d m u 




(i + yl K ) 


l±K 


m a m u 




(i-v 12 vl-yi K ) 


e7r, erj 


m d m u 




o-vuv^-yf) 


A i7r > M 7 ? 


m s m u 







Taking all these factors into account one can arrive at the following rough 
hierarchy of branching ratios. 

BR(vK) > BR(Dtt) > BR(IK) > BR(ln) (18) 

The most dominant decay mode of the proton normally is DK + . This 
pattern can be modified in some situations because of the interference of the 
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contributions from the third generation vs second generation making Dtt + . 
the most dominant decay mode. However, aside from this special situation 
which may correspond to fine tuning, the most dominant decay will be i?K + . 
We discuss now the details of this decay mode in minimaLsupergravity. The 
p — > i>iK + decay width for the neutrino type Vi is given bjtll 

T(j>^9 i K+) = (-§^) 2 \A\ 2 \B i \C (19) 
M H3 

Here Mh 3 is the Higgs triplet mass and (3 P is the matrix element between the 
proton and the vacuum state of the 3 quark operator so that 

(3 p Ul=e abc e a0 <OK L 4 L u1 L \p> (20) 

where U2 is the proton spinor. The most-reliable evaluation of (3 P comes from 
lattice gauge calculations and is given b\E!l 

I3 p = (5.6 ± 0.5) x lO- 3 GeV 3 (21) 

The other factors that appear in Eq.(19) have the following meaning: A con- 
tains the quark mass and CKM factors, Bi are the functions that describe 
the dressing loop diagrams, and C contains chiral Lagrangian factors which 
convert a the Lagrangian involving quark fields to the effective Lagraingian 
involving mesons and baryons. Individually these functions are given by 

2 

A = 5 ^m a m c t£ 1 VMiA S f (22) 

ZM W 

where Vij are the CKM factors, and Al and Ar are the long distance and the 
short distance renormalization group suppression factors as one evolves the 
operators from the GUT scale down to the electro-weak scale. Bi are given by 

_ 1 mfV^ rp p m t V 31 V 32 

u i— ■ 0f j — r-r 2 ±*2iH 7} T/ ^3B 3i \ (26) 

sm2(i rrisV^ w c V21 V22 

where the first term in the bracket is the contribution from the second gener- 
ation and the second term is the contribution from the third generation. The 
functions Bji are the loop intergrals defined by 

Bji = F(v,i, dj, W) + {dj -► e s ) (24) 

where 
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F(v,i, dj, W) — [Ecos^i_sin^ + f{ui, dj, W\) + cos^ + sin , y_f{u i , dj, W\)] 

1 8 l3 m"sin25 ui . . ~, ~ - t t ; n 

~o 7^T7 — [£/s«n7_stn7+/(iiii,dj,VVi) - cos^cos^+fiun, dj, W 2 ) 

-(Oil -> «i2)] (25) 

and / is given by 

f(Ui,dj, Wk) = sin 2 S ui f(ua,d j ,Wk) + cos 2 S u J(ui2,dj 1 W k ) (26) 

with 



/(a, 6, c) = 2 mc 2 [ 2 m " 2 fn(^4) - (m„ - m c )] (27) 



7 ± = /3+ ± P. 



and 



V2^± = M w (sin(3 ± cos/5) (29) 

sin2d U 3 = 5 ^ (30) 

m~ — mf 
ti t 2 

E = 1 , sm2/3 > fj,m 2 /M^ 
= -1, sin2/3 < iJ,rh 2 /M^ (31) 



Finally C is given by 



C=^[(l + ^±^)(1-^)? (32) 
327r/^ tub m z N 

where /„■, D, F, .. etc are the chiral Lagrangian factors and have the numerical 
values: U = 139 MeV,D=0.76,F=0.48,TOjv=938 MeV, m K =495 MeV, and 
m B =1154. 

The proton can also decay into vector bosons and these modes consist of 

eK* , fiK* ,ep, np,euj, fiuj (33) 
A similar analysis can be done for these mode@. 



Wc discuss now details of the analysis in minimal supergravity unification. 
In the analysis one looks for the maximum lifetime of the proton as we span 
the parameter space of the minimal model within the naturalness constraints 
which we take to be 

too < ITeV, mi < ITev, tan/3 < 25 (34) 

In Fig. 1 we display the results of the maximum lifetime for the mode p — > vK + 
as a function of the gluino mass. The analysis indicates that if Super^K and 
Icarus can reach the maximum sensitivity of 2 x 10 32 y for this modfflcJ then 
most of the parameter space within the naturalness limits indicated will be 
exhausted. The analysis including the dark matter constraints was also car- 
ried out. In this analysis we impose the very conservative constraint that the 
mass density of the relic neutralinos in the universe not exceed the critical 
relic density needed to close the universe. The results are also displayed in 
Fig 1. Here one finds that this constraint limits the maximum gluino mass to 
lie below 500 Ge\o. Most of this gluino mass range can_be_explored at the 
upgraded Tevatron with an integrated luminosity of abouOc3 25/6 -1 . 

5. Nucleon Stability in Extensions of the Minimal Model 



In Sec. 4 we discussed the situation regarding proton decay in the minimal 
SU(5) supergravity model. We discuss now the situations in some extensions 
of the minimal model. There are various kinds of extensions of the minimal 
SU(5) model. These include extension to larger groups, inclusion of textures, 
flipped models, string models, etc. We will consider here mostly the first two. 
In SO(10) one has a large tan/3, i.e., tan/3 ~ 50 to achieve b — t — r unification 
and compatibity with the observed mass ratios for b, t and r. However, a large 
tan/3 tends to destabilize the protorE3. This can be seen from the fact that the 
current lower experimental limits for the proton lifetime for this decay mode 
impose the following constraint on the effective mass scale Mpp> = (M _1 )n 

M PD > tanf3(0.57 x 10 16 )GeV (35) 

?w) is plotted in Fig. 2 taken from refEH. One finds that for 
values of Mpp, of O(10 17 ) GeV, the disagreement oLthe theoretical value of 
a s in minimal SO(10) with experiment is about 6c£3. Thus one needs_iarge 
threshold corrections at the GUT scale to achieve agreement with dat£j. 

Next we turn to another type of extension of the minimal SU(5) model 
and it involves inclusion of textures to generate the correct quark-lepton mass 



which for tan/3 = 50 gives Mpo ~ 2.5 x 10 17 GeV. However, this large scale 
tends to upset the unification of the gauge coupling constantOo. The analysis 
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300.0 500.0 700.0 900.0 

Gluino mass (GeV) 



Figure 1: Plot of the maximum p — » PK + lifetime under the constraint mo < 1 TeV in 
minimal supergravity model. The solid curve is for the case when no relic density constraint 
is imposed. The dashed curve is with the dark matter constraint Qh 2 < 1. The dashed 
horizontal line is the current experimental limit and the solid horizontal line is the limit 
expected at Super-K and Icarus experiments. 
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0.2308 


0.231 


0.2312 0.2314 0.2316 0.2318 
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Figure 2: Contour plot of in sin 2 (9w) - a^(Mz) plane. The dashed line corresponds 

to the lower bound on M PD in minimal SU(5) of 1.2 X 10 16 GeV. The dot-dashed line 
corresponds to the lower limit on Mpd in minimal SO(10) of 2.7 X 10 17 GeV. Currentlyprtie 
measurements are sin 2 (0 w ) = 0.2313 ± 0.0003 and a 3 (M z ) = 0.119 ± 0.0003 (from refEJJ). 
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hierarchicsE! One procedure to generate the textures is to include Planck 
scale corrections in the superpotentiaEiE3, where the Planck scale corrections 
involve expansions in the ratio Yl/Mpi, and £ is the adjoint scalar field. After 
spontaneous breaking of the GUT group £ develops a VeV and one gets a 
hierarchy in the ratio M/Mpi which allows one to generate the textures in the 
Higgs doublet sector defined by 

W d = H 1 lA E e c + H 1 d c A D q + H 2 u c A u q (36) 

where A E , A D and A u are the textures matrices. Once the textures in the 
Higgs doublet sector are fixed one can compute the textures in the Higgs triplet 
sector defined by 

Wt = H 1 lB E q + H 2 u c B u e c + t abc {Hid c b B D u c c + H^u c b C u d c ) (37) 

where B E , B u , B D , and C u are the triplet textures. 

For the case when one adds the most general Planck scale interaction one 
finds that fixing the textures in the Higgs doublet sector does not uniquely fix 
the textures in the Higgs triplet secto£3. One needs a dynamical principle to do 
so. One possibility suggested is to extend supergravity models ta include not 
only the usual visible and hidden sectors, but also an exotic secto£3. The exotic 
sector contains new fields which transform non-trivially under the GUT group 
and couple to the fields in the hidden sector and the adjoint scalars in the visible 
sector. After spontaneous supersymmetry breaking the exotic fields develop 
Planck scale masses because of their couplings with the hidden sector fields and 
integration over the exotic fields leads to the desired Planck scale corrections. 
For the choice of a minimal set of exotic fields one finds that the Planck scale 
corrections are uniquely determined and one finds that correspondingly the 
textures in the Higgs triplet sector are uniquely determined. 

We give now some specifics. The simplest examples of textures are the 
Georgi-Jarlskog matrices given by 

/ F \ / Fe l(t > \ / C \ 

A E = I F -3E , A D = Fe'^ E \ ,A U = \C B J 
\0 D ) \ D J V B A) 

(38) 

If one uses the most general Planck scale interaction then fixing the texture 
in the Higgs doublet sector to be the Georgi-Jarlskog does not fix the textures 
uniquely in the Higgs triplet sector. However, for the case when one uses the 
exotic sector hypothesis with the minimal set of exotic fields one finds that the 
textures in the Higgs sector are uniquely determined and are given bjEHl 
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/ aF \ / —§fF 

B E =\a*F fE ),B D = [ ^F -§£ j (39) 
V \d) \ -\D t 

/ |C \ 
B^ 7 = |C -Ifl (40) 

\° -\ B A J 

where a=(— i= + e 1 ^) and C* 7 = _B y . Estimates including the textures show 
modest (i.e. 0(1)) modifications in the p decay lifetimes. Further, the various 
p decay branching ratios are affected differentially and thus measurement of 
the branching ratios will shed light on the textures and on the nature of physics 
at the GUT scale. 

Finally, we discuss briefly proton stability in string models. It is difficult 
to make general remarks here as there are a huge variety of string models 
and further the mechanism of supersymmetry breaking in string theory is not 
fully understood and thus the nature of soft SUSY breaking parameters is not 
known. However, there are various types of parametrizations that have been 
used in the literature. One such parametrization is that of no-scale models 
where one assumes that at the GUT scale one has mo = = Aq and radia- 
tive breaking of the electro- weak symmetry is driven by mi and the top mass 
rrit- In this case it is shown that the minimal SU(5) would lead to a proton 



decay via dimension five operators whic 
experimental bounds for the vK + mod 



n would be in violation of the existing 
This difficulty is eliminated in the 



flipped no-scale models where the DK + mode is highly suppressed. Of course, 
as mentioned above there is a huge array of string models and each model must 
be individually examined to test p stability. Thus an analysis of p stability in 
three generation Calabi-Yau models can be found in reO. Similarly an analy- 
sis of p stability in superstring derived standard like models is given in ref. l3. 
7. Conclusion 



In this review we have given a brief summary of the current status of proton 
instability in SUSY, SUGRA and string unified models. One of the important 
conclusions that emerges is that if the Super-K and Icarus experiments can 
reach the expected sensitivity of 2 x 10 34 y for the vK + mode, then one can 
probe a majority of the parameter space of the minimal SUGRA model within 
the naturalness constraint of mo < 1 TeV, m~ g < 1 TeV, and tanj3 < 20. If one 
includes the additional constraint that the neutralino relic density not exceed 
the critical matter density needed to close the universe then the gluino mass 
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must lie below 500 GeV to satisfy the current lower limit on the vK + mode. 
These results will be tested in the near future in p decay experiments as well 
by experiments at the Tevatron, LEP2 and the LHC. 
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